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Abstract 

The purpose of this paper is to study quasinormal modes (QNM) of a regular black 
hole with a cosmological constant due to scalar perturbations. A detailed study of the 
QNM frequencies for the massless scalar held were done by varying the parameters 
of the theory such as the mass, magnetic charge, cosmological constant, and the 
spherical harmonic index. We have employed the sixth order WKB approximation to 
compute the QNM frequencies. We have also proved analytically that the I = 0 mode 
for the massless held reach a constant value at late times. We have approximated the 
near-extreme regular-de Sitter black hole potential with the Pdschl-Teller potential 
to obtain exact frequencies. The null geodesics of the regular-de Sitter black hole is 
employed to describe the QNM frequencies at the eikonal limit {I >> 1). 
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1 Introduction 

It is a well known fact that our universe is expanding at an accelerated rate. There 
are many observational evidence for this [I] [2] |3] |1] [5] . One of the simplest candidates 
to describe the mysterious dark energy causing the acceleration is the cosmological 
constant A introduced by Einstein to General Relativity to describe a static universe. 
The space-time with a positive cosmological constant is known as the de-Sitter space- 
time. Hence, our universe may be described by de Sitter geometry. There are several 
reasons to study de Sitter spaces. In high energy physics, there is a holographic 
duality relating quantum gravity on the de Sitter space to conformal held theory on 
a sphere of one dimension lower 0 0 known as dS/CFT correspondence. Examples 
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along those lines include, Kerr-de Sitter space-time in various dimensions studied 
by Dehghani |8] and holographic description of the three dimensional Kerr-de Sitter 
space with gravitational Chern-Simons term studied by Park [9]. One other reason for 
de Sitter space-times to gain lot of attention is the need to understand it in from the 
point of string theory. String theory is one of the popular candidates for a theory of 
quantum gravity. Hence, a fully satisfactory de Sitter solution to string theory has to 
be found. Towards this goal, classical solutions in Type IIB string theory have been 
studied by Dasguptha et.al [10]. Another work studying classical de Sitter solutions 
in string theory has been published by Danielsson et.al. m- 

Usually, a black hole is a space-time with a horizon and a singularity covered by 
the horizon. Contrary to this notion, there are black holes where there is a horizon 
but without a singularity at r = 0. Such black holes are called “regular” black holes. 
The first regular black hole was constructed by Bardeen [T2| [T3| and interpreted as a 
magnetic monopole by Aydn-Beato and Garcza [H]. There are several other regular 
black hole solutions proposed by Aydn-Beato and Garcia ransmiiiia. In this 
paper, we are interested in studying one of the most interesting regular black holes 
proposed by Aydn-Beato and Garcia [15] and inteprited by Bronnikov [19]. This black 
hole is static, spherically symmetric and is described by two parameters, mass M and 
magnetic charge Q. For >> 1 and for small such black holes have similar 
properties as the well known Reissner-Nordstrom solution. This family of black hole 
solutions are referred to as ABGB black holes. Here we study the ABGB black hole 
in a de Sitter universe presented by Matyjasek et.al [SOj. Since the current universe 
may be a one with a positive cosmological constant, it is important to understand 
regular black hole embedded in a de Sitter universe. 

When a black hole is externally perturbed, the response can be described by three 
stages. The initial stage entirely depend on the initial data of the external perturba¬ 
tion. The intermediate stage correspond to damped oscillations with incomplete set 
of complex frequencies. These modes are called quasi-normal modes (QNM) and de¬ 
pend only on the macroscopic properties of the black hole such as the mass, spin and 
the charge. Measuring QNM frequencies using current gravitational wave detectors 
such as LISA, VIRGO and LIGO give a direct approach to detect black holes pi] . 
Other than for experimental reasons, there are also several theoretical aspects of black 
holes which are studied by using QNM frequencies. Due to the famous AdS/GFT 
correspondence, many works have been devoted to study QNM’s of black holes in 
anti-de Sitter space [22] [23] . Black hole quantum spectrum via QNM were presented 
by Gorda in [21]: the black hole model used by Gorda is somewhat similar to the 
semi-classical Bohr’s model for the hydrogen atom. Another reason to create inter¬ 
est on QNM’s was the conjecture by Hod [25] relating quantum properties of black 
holes and asymptotic values of frequencies of QNM’s. There are many works de¬ 
voted to compute asymptotic QNM frequencies due to this reason. Some recent work 
have addressed relation between QNM’s and hidden conformal symmetry of wave 
equations on the black hole back ground [2S] [2Z]- Other than black holes, naked 
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singularities also have been studied from the QNM point of view [28]. An excellent 
review on QNM’s is written by Konoplya and Zhidenko [29]. Given all the above, it 
is worthwhile to study QNM’s of regular black holes with nonlinear sources. 

The paper is organized as follows: in section 2 an introduction to the regular-de 
Sitter black hole is given. In section 3 the scalar held perturbation is introduced. 
In section 4, the WKB approach to compute QNM are introduced. Massless scalar 
perturbation with QNM frequencies are presented in section 5. In section 6, analysis 
of / = 0 mode of the massless scalar held is presented. In section 7, unstable null 
geodesics approach is used to compute QNM frequencies. In section 8, Pdschl-Teller 
method to compute QNM frequencies are presented. In section 9 the conclusion is 
given. 


2 Regular black hole with a cosmological constant 
(ABGB-de Sitter black hole) 

In this section we will present the basic characteristics of a regular black hole with 
a positive cosmological constant derived by Matyjasek et.al [20]. The action from 
which it was derived is, 


S = -^I^HR-2A)-C(F)] ( 1 ) 

Here R is the scalar curvature and A is the positive cosmological constant. C{R) is 
the Lagrangian for non-electrodynamics component given by. 


C{F) = F 


1 — tanh?' 



( 2 ) 


Here « = ^ and F = Ff^^F^’^ with = SJ The only component of 
Ffj^ty which is non-zero for this solution is F 23 = Fg^ given by. 


F 23 = QsinO 


Here, Q is the magnetic charge. Hence, 

rpA. 


( 3 ) 

( 4 ) 


Static spherically solutions for the action in eq.([T]) with the given magnetic held was 
derived in [20] to be. 


ds^ = -f(r)dt^ + 

fir) 


( 5 ) 
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where, 


2M 


Ar^ 


f{r) = 1 - 


\2Mr 


( 6 ) 


This particular solution was named Aydn-Beato-Garcza-Bronnikov-de Sitter (ABGB- 
dS) solution. 

The function /(r) is finite at r = 0. /(r) can have at most three distinct roots 
leading to three horizons: Vc is the cosmological horizon, rt is the event horizon and 
Tj is the inner horizon. The region considered for the perturbations will be between 
Tfe and Vc- f{r) also can have degenerate horizons: when they are referred to 

as cold black holes and when rf, = rc, they are considered as charged Nariai black 
holes [30] [31]. 

When Q = 0, the metric yeilds the Schwarzschild-de Sitter black hole and when 
A = 0, one obtain the ABGB black hole. Function /(r) can be expanded for small Q 
and for large r as. 


2M ^ ^ q6 

r 3 12M^r‘^ 


(7) 


Hence, it is clear that for large r and for small Q, the line-element resembles the well 
known Reissner-Nordstrom-de Sitter black hole (RNdS) with the line element eq.f[Sl) 
with the function /(r). 


f{r)iiNdS = 


2M 

r 3 


( 8 ) 


The horizon structure of the ABGB-de Sitter black hole has been analyzed thor¬ 
oughly in [20]. The regularized stress-energy tensor of the scalar, spinor and vector 
fields inside the degenerate ABGB-de Sitter black hole solution were studied by Maty- 
jasek et.al PI. 


3 Scalar field perturbations 

Having introduced the geometry of the black hole space-time, in this section we will 
present perturbations by a minimally coupled scalar field around this geometry. The 
equations of motion for such a held is given by the Klein-Gordon equation, 

(V^ - = 0 (9) 

The scalar held cp can be decomposed into its partial modes as 

r]{r, 9, (j),t) = Y^ e~^‘^%^rn{9, (p) (10) 

, r 

l,m 

Here, Yi^rn{9, (p) are spherical harmonics, and, I and m are the angular and the mag¬ 
netic quantum numbers, oo is the frequency of the oscillations of the held. With 
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the above decomposition, the radial component of the Klein-Gordon equation can be 
simplihed to be, 

- Vscalarin)) R{n) = 0 ( 11 ) 

Here r* is the tortoise coordinates which is given by, 

dr. = ^ ( 12 ) 

f{r) 

and Vscaiar is given by, 

= /(r) + ^ + y) (13) 


The effective potential approaches zero at r = ad r = Tc. It has a peak in between 
as demonstrated in Fig.([T]). 

To get approximations for r., one can expand /(r) around r = in a Taylor 
series as, 

f{r) ^ f{rb,c){r -rb,c) (14) 

Since f'{rc) < 0, r. can be derived to be, 

r* - ~ |/^(r )\ ^^d{rc - r) (15) 


Hence, when r —>• Vc, r. —)■ oo. Similarly, since f'{rb) > 0, r. can be derived to be, 

1 


\f'i.rb)\ 


Log{r - Tb) 


(16) 


Hence, when r —)■ rb, r. —oo. Hence the effective potential Vgcaiar —t 0 when 
r* ±cx). 


4 Computation of QNM frequencies of the scalar 
perturbations by WKB approach 

QNM for a black hole perturbed by a scalar held is given by the solutions of the 
wave eq. flllD . To hnd the QNM solutions one has to impose boundary conditions; 
it is purely ingoing held at the even horizon and purely outgoing helds at the 
cosmological horizon r^. The frequencies corresponding to the QNM are given by 
Lo = Ljfi + iui; here ur is the oscillating component and w/ is the damping component 
of the frequency. QNM for eq. flTTl) with the above mentioned boundary conditions 
can be represented as, 

R{r.) exp{iujr.); r. —oo{r Vb) (17) 
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-R(r*) —)■ exp{—iur^); r* +co(r Tc) (18) 

Since the effective potential has a peak and has the form of a potential barrier, to 

hnd the QNM frequencies, one can use the WKB approach. WKB approach to find 
QNM frequencies of black holes was hrst developed by Iyer and Will for third order 
[32] : it was later extended to sixth order by Konoplya [33]. The sixth order WKB 
approximation is employed to find QNM frequencies in several papers including, [13] 

m ra- 

In the sixth order WKB approach, QNM frequencies are given by the expression, 

21/"(rmaa;) ^ 2^ + ^(19) 

Here, Vmax is where V{r) reach a maximum and V''{r) is the second derivative of the 
potential. Expressions for Tj can be found in [33] . 

5 Massless scalar field perturbation 

First, we will focus on the massless scalar held perturbations where /x = 0. In this case, 
the effective potential Vscaiar depend on the parameters, M,QJ and A. In Fig.([T]), 
the potential is plotted by varying M. When the mass M increases, the height of the 
potential decreases. In Fig.(|2]), Vgcaiar is plotted by increasing the charge Q. When Q 
is increased, the height of the potential increased. When the cosmological constant A 
is increased, the height of the potential decreases as given in Fig. ([3]). In Fig(|lj), the 
potential is plotted by varying I where I is positive. When I increases, the potential 
increases. In Fig([5]), the potential is plotted for I = 1 and / = 0. It is interesting to 
notice that the potential for I = 0 has a local minimum between the two horizons, 
and Tc which is different from the behavior of the potential for other values of /. 


V(scalar) 



r 


Figure 1: The figure shows Vscaiarij") vs r. Here Q = 0.086, A = 2.5 and I = 2 


6 





14 
12 
10 

V(scalar) ^ 

6 

4 
2 
0 

0.2 0.3 0.4 0.5 0.6 



Figure 2: The figure shows Vgcaiarir) vs r. Here M = 0.1222, A = 2.5 and / = 2 
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Figure 3: The figure shows Vscaiar{f’) vs r. Here Q = 0.1055, M = 0.129 and / = 2 
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V(scalar) 



Figure 4: The figure shows Vscaiarij") vs r. Here Q = 0.1055, A = 1.835 and M = 0.129 


V(scalar) 



Figure 5: The figure shows Vscaiar{r) vs r for / = 1 and I = 0. Here Q = 0.1055, A = 
1.835 and M = 0.129 


5.1 Quasinormal modes of massless scalar perturbations for 
I > 0 modes 

In this section, we will discuss the properties of QNM frequencies for massless scalar 
fields for / > 0 modes. We want to mention that for all u computed ojj was negative 
leading to the conclusion that the black hole is stable under scalar perturbations for 
/ > 0. When plotting, we will only use the magnitude of ui in the rest of the paper. 

First, we computed QNM frequencies by varying charge Q] ojr vs Q and ui vs 
Q are given in Fig.(|6]) (for I = 2) and Fig. ([8]) (for / = 1). When Q increases, ujr 
increase. On the other hand when Q increase, coj increase to a maximum and then 
start to fall off at large Q. Hence, there is a maximum value of Q where the black 
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hole has the most stability. We want to note that when Q become large, the black 
hole is approaching the degenerate limit where r^. Also, we have compnted the 
“Qnality factor” for the held which is dehned as, 

(jJ p 

Quality Factor = ——r (20) 

2|w/| 

If the qnality factor is large, then the held is considered a better oscillator. The 
qnality factor for the scalar held is plotted in Fig. ([7]) (for I = 2) and FigdH]) (for 
/ = 1). The larger the charge, the larger the qnality factor. Hence, the scalar helds 
are better oscillators for large charge. 

We compnted the radins of the black hole horizon for the same valnes of Q and 
plotted u vs Tb in Fig. flT0|) . First we want to mention that Vb decreases when Q is 
increased. Hence it is no snrprise to see that ojr decreases when rb increases. On the 
other hand, the behavior of ojj vs Vb is similar to the behavior of uj vs Q: there is a 
particnlar valne of Vb where ojj reaches its maximnm which is the most stable radius 
of horizon where stability is concerned. 

For the same Q values, we computed the Hawking temperature of the black hole; 
the plot of uj vs Tb is given in Fig. fllip . In many black holes, uj seems to have a 
linear relation with the temperature [23] • In the regular black hole, there seems to 
be a strange relation between uj and T; it is possible to have two values of ojj for the 
same temperature. 

Next, we computed oj by varying A. The plots for ur vs A and uj vs A are plotted 
in Fig. ffT2|) . Both ur and uj decreases as A increases. Hence, a smaller A is favored 
for a held decaying faster around the black hole. The quality factor as a function of 
A is plotted in Fig. ffT31) . One can see that the held is a better oscillator for small A. 

QNM frequencies were also studied by varying 1; ojr vs I is plotted in Fig. ffTT|) and 
uji vs I is plotted in Fig. ffTSj) . When observed the behavior of cur, it is clear that it 
varied linearly with I for both n = 0 and n = 1. When observed coj, it is clear that it 
decreases to reach a stable value when I increases, for both n = 0 and n = 1. Hence, 
the held decays faster for small I values. 

We have also plotted cu vs n in Fig. ffTbj) . Since the accuracy of the WKB approach 
is greater for I > n, we have chosen I = 10 for this computation. In this case, ur 
decreases with n which is opposite to what we observe in a simple oscillating system 
with standing waves such as a tight string. On the other hand, uj increases linearly 
with n. 
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Figure 6: The figure shows u vs Q. Here M = 0.9, A = 0.076 and I = 2 


Quality factor 



Figure 7: The figure shows Quality factor vs Q. Here M = 0.9, A = 0.076 and I = 2 
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Figure 8: The figure shows u vs Q. Here M = 0.9, A = 0.076 and I = 1 



Figure 9: The figure shows Quality factor vs Q. Here M = 0.9, A = 0.076 and I = 1 
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Figure 10: The figure shows u vs r?,. Here, M = 0.9, A = 0.076 and I = 2 
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Figure 11: The figure shows oji vs Tb. Here, M = 0.9, A = 0.076 and I = 2 




Figure 12: The figure shows u vs A. Here Q = 0.25, M = 0.278 and I = 2 
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Figure 13: The figure shows Quality factor vs A. Here Q = 0.25, M = 0.278 and 
I = 2 



Figure 14: The figure shows ojr vs 1. Here, Q = 0.25, M = 0.278 and A = 1 
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Figure 15: The figure shows ui vs 1. Here, Q = 0.25, M = 0.278 and A = 1 




n n 

Figure 16: The figure shows u vs n. Here Q = 0.25, M = 0.278 ,1 = 10 and A = 1 


6 Analysis of / = 0 mode of massless scalar field 

As given in Fig.(|5]), potential for / = 0 modes is rather different from all other potential 
with / > 0. The difference is attributed to the fact that the potential is negative in a 
range between the black hole horizon and the cosmological horizon. In fact, similar 
behavior persists for I = 0 potential for both Schwarzschild-de Sitter black hole and 
the Riessner-Nordstrom-de Sitter black hole as shown in Fig flTTj) . 
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Figure 17: V{r) vs r for I = 0 mode for Reissner-Nordstrom-de Sitter and 
Schwarzschild-de Sitter black hole. Here, M = 0.1082, A = 2.5 and Q = 0.086 



WKB approximation that was used in section(5.1) to compute the QNM frequen¬ 
cies of / > 0 modes are not suitable to compute the frequencies of / = 0 modes since 
there is an extra local minimum of the potential. Hence we will avoid computing the 
QNM frequencies for this case. Instead, we would be interested in understanding if 
the black hole is stable under the I = 0 mode perturbations. In fact, when the po¬ 
tential is negative in some region, it is possible to have growing QNM modes leading 
to instability of the system under such perturbations. However, it turns out that this 
is not the case for all cases with negative potentials as discussed in [36]. What we 
would present here is to the fact that the held of / = 0 mode in fact settle in to a 
constant value at late times. Since one of the reasons to study QNM frequencies is 
to understand the stability of the black hole under such perturbations, such a study 
is appropriate. If one needs to compute the QNM frequencies, and see how the held 
evolves, then a time-domain integration method used by Gundlach et.al. in [3H] can 
be used. Instead, we will present an analytical observation of how the held behaves 
at late time. 

I = 0 mode does approach a constant value for the Schwarzschild-de Sitter black 
hole as presented by Brady et.al [37|. We will use similar approach to show the 
behavior of the held around the regular-de Sitter black hole. 

Here, we will study the scattering of the held between and Tc for I = 0 mode. 
First we will introduce null coordinates u and v as, 

u = t — r^] v = t + r^ (21) 

Here, r* is the tortoise coordinate introduced in eq. llT^ . Then, the line element in 
eq.([5|) simplihes to 

ds^ = —f{r)dudv + r'^dVt^ (22) 
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According to the null coordinates, the future cosmological horizon, r = Tc is located 
at u = oo and the future black hole horizon, rj, is located at u = oo. Now, in terms 
of the null coordinates, the scalar wave equation X/'^rj = 0 for a massless scalar held 
simplihes to. 


(9^\l'(u, v) 
dudv 



(r)'h(u, v) 


(23) 


where, 

r]{u, V, e,(j)) = Y^ 0) (24) 

l,m 


and 


4^=0 = 


r 


(25) 


Notice that for / = 0, the potential simplihes. 

When / = 0, the solution to the eq. fl2^ can be expressed as a series depending on 
two arbitrary functions G{u) and H{v) as [37][38] . 


^ = A (G(m) + H{v)) + ^ Bp{r)[G^-P-^\u) + H (26) 

p=0 


Here, the negative super indices on G and H refer to integration with respect to u 
and V respectively: for example, G“^(u) = / G{u)du. The coefficient A can be set 
equal to 1 without loss of generality. 

When eq. fl26ll is substituted to eq. fl2^ . one can compare the coefficients to come 
up with a recurrence relations for Bp{r) as, 

= h (_B^ + ,b;) + Lb; (27) 

Here, 

^0 = - / ^dr (28) 

To facilitate the computation, we will assume the charge Q is small and do an expan¬ 
sion for the function /(r) for the regular de-Sitter black hole as. 


fi.r) 


2M 

r 



12M2r4 


(29) 


Note that the qualitative behavior of the black hole outside and inside are not 
much diherent in the expansion. Now, one can integrate eq. fl28|) to obtain Bq as. 


M rA 

° 30M2r5 3r3 ^ 2r2 ^ 3 

From the recurrence relation in eq. flTTjl . Bi can be obtained to be. 


(30) 


Bi = 


11 


E 


(Xfi 

j-n 


(31) 
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where an depends on M, Q and A. Other Bp values with p > 1 also behaves similarly 
with higher powers of ^ terms. For large r, Bp will vanish except for Bq. This is one 
of the major difference compared to an asymptotically flat black hole where all Bp 
values will vanish for large r. 

Now, we consider an initial burst of radiation on n = 0. The radiation is conhned 
between u = 0 and u = Ui and is given by 0) = G{u). Also, \h(0, v) = 0. This is 
demonstrated in Fig. ffTSjl . We want to study the evolution of this held in the diamond 
shape region given in Fig. fllSI) by thick lines. 



Figure 18: The region in space-time outside the black hole studied for the evolution 
of an initial burst conhned between u = 0 and u = Ui and on n = 0. The initial bust 
is such that, T(tt, 0) = G{u) where G{u) 7 ^ 0, for 0 < ti < Mi only. 


Following the work by Gundlach et.al [38] and Brady et.al. [3^, the evolution 
of the initial burst can be broken into two parts: (i) the evolution of the burst 
T(m, 0) = G{u) which is non-zero between 0 < u < ui and \k(0, v) = 0. (ii) subsequent 
evolution in the region for u > Ui. 
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6.1 The evolution of the burst = G{u) which is non¬ 

zero between 0 < u < mi and ^'(O, v) = 0 

Since the initial burst is independent of v and G(mi) = 0, the held a.t u = Ui is given 
by, 

oo 

= Y.Bp{r)G{uY-P-^\u^) (32) 

p=0 

For sufficiently small A, Vc is large and at r = the only non-zero S is i?o ~ 
Hence, at the cosmological horizon Vc {v = oo), 

Ar 

'^{ui,oo) ^ BoG~^{u) = — - G{u)du + 0{A) (33) 

3 Jo 


6.2 The evolution of the field in u > ui region 

Now, we can examine the evolution of the held in the region u > uihy taking eq.(j32]) 
as the initial data. Following the work in |37] and [3H], the held equation given in 
eq. fl2^ can be solved near the cosmological horizon as, 

^{u,v) ^ [G{u) + H{v)] + Bo{r)[G-\u) - H-\v)] (34) 


Dehning G ^{u) = g{u) and H ^(u) = —h{v) where g{u) and h{v) are arbitrary 
functions. Approximating Bq to be ^ since r is large, one can rewrite eq. fl3^ as. 


4'(w,n) —[g{u) + h{v)] + 

Since closer to the cosmological horizon r^, r* 
also the fact that r* = one can write r as. 


dg{u) dh{y) 
du dv 

'7W)^°s(rc- 


(35) 


r) from eq. lfTSj) . and. 


r 


j. _ ^-KcV^KcU 


(36) 


Here v = 


df'(r) 

dr 


is the surface gravity at Tc. Hence closer to the cosmological 

c 

horizon (when v —>■ oo), the right hand side of eq. fl32p can be approximated with a 
polynomial of Hence by comparing eq. (l5^ and eq. (l5H]) . one can conclude that 

the function h{v) is also a polynomial of given as. 


h{v) = E 

71=0 


(37) 


When n ^ oo, one can compute Hq from eq. fl33|) and eq. (l35D . as, 

3 

ho = ^'h(Mi,oo) 


(38) 
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which is non-zero. The nature of the function g{u) depends on the potential every¬ 
where. However, since g{u) also evolves from eq (l3^ . it seems reasonable to predict 
that it is of the form (when u —)■ cxd) 


OO 

g{u) = Y.gne-^^^^ (39) 

Hence, at late times, the field T —» ho+go and the scalar field g = ^ = approach 
a constant value at late times (when u and v both approach oo). This is only possible 
only when ho + go ^ 0. When I = 0 mode, this is the only static solution which is 
regular at both horizons and r^. Hence, the field decay into a constant value late 
times only for this particular mode without decaying completely. Therefore, I = 0 
mode of the massless scalar field is prone for instability. 

7 Quasi normal modes of massless fields from null 
geodesics 

The black holes considered in this paper have unstable null geodesics. Cardoso et.al. 
|39] . showed that in the eikonal limit (/ >> 1), QNM of a black hole can be determined 
by the parameters of null geodesics. An earlier work done along those lines is given 
by Mashhon [ID] 

An expansion method to compute QNM based on the link between the null 
geodesics and QNM freuencies has been applied to Kerr black hole in a paper by 
Dolan [H]. The method of null geodesics has been applied to study the near-extreme 
Kerr black hole by Hod |12] and black holes in anti-de Sitter space by Morgan et.al 
in [13]. 

The equation of motion for massless particles around a black hole is given by [H] 

+ Vnuii = (40) 

with, 

VnuU = /(r) (41) 

Here, L is the angular momentum of the massless particle. For bot r = ru and r = Tc, 
Vnuii = 0. In the Fig. ([T9]), Vnuii is plotted for various values of the magnetic charge 
Q. One can observe that the height of the potential is higher for higher values of Q. 
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Figure 19: The figure shows Vnuii{'>') vs r. Here M = 0.19, L = 1 and A = 2 


Now, when the effective potential for the massless scalar field Vscaiar in eq.(17) is 
studied, one can see that in the eikonal limit (/ —)■ cxd), 

Vscaiar ~ ^ (42) 

Hence, one can conclude that the maximum of Vscaiar occurs at r = Vmax given by, 

2/'(r’maa:) l^maxf ij" max) 0 (43) 

On the other hand, the maximum of Vnuii occurs at V^^n = 0 leading to, 

‘^firnull) - rnull f'{rnull) = 0 (44) 

Hence the maximum of Vscaiar and the location of the maximum of the null geodesics 
coincides at r^ax = Tnuii- Vnuii{r)' = 0 leads to the equation, 

%Sech^ {+ 4 f-3M + r + 3MTanh {^ = 0 (45) 

r5 \2Mr) \2Mr)) ^ ' 

Solutions to the above equation gives Vnuii values. Since an analytical expression is 
not possible, we have calculated Vnuii varying charge Q numerically and plotted in 
Fig. fl20l) . One can observe that Vnuii gets smaller as Q increases. 
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Figure 20: The figure shows VnuU vs Q. Here M = 0.1715, L = 1 and A = 0.0973 


According to [39] the QNM frequencies in the eikonal limit is given by, 

^^QNM = +-)\XLy\ (46) 

Here, I is the angular momentum of the perturbation and n is the overtone number, 
and Ha is the coordinate angular velocity given as, 


Ha = 


i 


(47) 


The parameter X^y appearing in eq. fldOp is the Lyapunov exponent which is interpreted 
as the decay rate of the unstable circular null geodesics. The derivation of the above 
results is explained in Cardoso et. al[39]. For the ABGB-dS black hole, Qa and Xiy 
are given as, 


i{rnuii) \ rl^ii 


(48) 


A 


Ly — 




2t(r„a«)^ 


null f null) 


2L2 


(49) 


By using Ha and the QNM frequencies are calculated and plotted in Fig. ([21]) 
and Fig.(|22|). When compared with u values computed with WKB approach given 
in Fig(|6|), one can conclude that they are similar in behavior. 
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Figure 21: The figure shows Xiy vs Q. Here M = 0.1715, L = 1 and A = 0.0973 



Figure 22: The figure shows vs Q. Here M = 0.1715, L = 1 and A = 0.0973 


8 Poschl-Teller approximation for the near-extreme 

ABGB-de Sitter black hole 

First we like to discuss the extreme ABGB-dS black hole when Vf, = Vc- For the black 
hole to be extreme, two conditions has to be satisfied: 

/(r) = 0; f{r) = 0 (50) 

Due to the complex nature of the function /(r), an exact value for for the extreme 
black hole is not possible. Let us name the value of Vf, at this special case as a. It 
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is clear that fie) < 0 due to the nature of the function f{r). Such extreme black 
holes where Vf, = are called Nariai black holes |31]|45][l6] . 

In the near extreme case, rt and Tc are very close. Hence /(r) can be expanded 
in a Taylor series as IMI as, 

f{r) ^ ^-^{r - rb){r - r^) (51) 

Hence the tortoise coordinate dehned as r* = / can be integrated as. 


dr 


r* = 


log 


fir) /"(cT)(rc - Tb) \r-rb 


(52) 


From the above derivation, it is clear that when r —)■ r^, r* —)■ —cxd and r —)■ Tc, 
r* —)■ oo. From the relation in eq. flS^ can be solved for r as. 


rb + 

1 + eP'^* 


(53) 


Here, 

P = - n) (54) 

Bn substituting r in eq. flS^ . the function /(r) for the near extreme ABGB-de Sitter 
black hole can be written as. 


fir) 


4{Cosh{^))^ 


Here, 





Now the effective potential for the scalar held can be written as. 


(55) 


(56) 




n 

Cosh^lf) 


(57) 


where 


Go = 


m+^) 

4(7^ 



(58) 


In deriving Vq, we have assumed r k, a. Also since f{r)' k, 0, only the hrst term in 
the potential dominates. With the new dehnition of the potential, the wave equation 
for the scalar held perturbations simplihes to. 


d?VL{r^) 
drl ^ 



Vo \ 
Cosh\^)) 


H(r*) = 0 


(59) 
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Cosf'^°( ) ' known Poschl-Teller potential and Ferrari and Maslihoon, in a 

well known paper m demonstrated that u can be computed exactly as, 


p- 


■P, 


oj = \ Vq - i—(n -\— 

V 16 2^ 2 


(60) 


By observing the expression for u given above, some hand waving arguments of the 
behavior of u can be made. When p is large, and / is small, Vq ~ Hence, 

ujji ^ + constant and ujr will depend linearly on On the other hand, uji does 

not depend on ^ for extreme black holes. When the spherical harmonic index I is 


large, W 


40-2 ■ 


Hence ur 




2a 


f constant. Hence ur depends on I linearly for large 


I which is clear from Fig flTT|) . For large I, oji is independent of I and is clear from 
Fig flT^ . For large n, ui becomes large and ur is independent of n. This behavior is 
represented in Fig. ffT^ . 


9 Conclusion 

Our main goal in this paper has been to study QNM and stability of a regular black 
hole in de Sitter universe under scalar perturbations. The regular black hole does not 
have a singularity at r = 0. It has three horizons: cosmological horizon Tc, black hole 
event horizon and inner horizon r*. The structure of the regular-de Sitter black 
hole is very similar to the Reissner-Nordstrom-de Sitter black hole. 

We have employed sixth order WKB approximation to compute QNM frequencies. 
The parameters of the theory, the mass M, magnetic charge Q, cosmological constant 
A and the spherical index I and the mode number n were changed to see how QNM 
frequencies depend on them. When Q is increased, ur increases. On the other hand, 
when Q is increased, uj increases to a maximum and then decreases for large Q. 
Hence, there is a particular value of Q where the black hole have the most stability. 
We have observed the behavior of the quality factor for the black hole as a function 
of the charge Q. It is noticed that the larger the charge, larger the quality factor is. 

Next, we have plotted a; as a factor of the black hole rf,. When increases, 
Ur decreases and uj increases to a maximum and then decreases. Hence there is 
a particular size of the black hole where the stability is maximized. We have also 
computed the Hawking temperature to observe the relation with ur The relation is 
not linear as in other black holes. 

We have studied the behavior of u with the cosmological constant A. When A 
is increased, both ur and ui decreases. Hence, a smaller A is favored for a stable 
black hole. The behavior of u with respect to I, the spherical harmonic index is 
similar to the behavior of other black holes, ur has a linear relation with I and uj is 
independent of I for large 1. We have also plotted u vs n, the mode number. Here, 
Ur decreases with n and uj increases linearly with n. 
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The / = 0 mode for the massless field presents a potential with a negative region 
between Vh and Vc- WKB approach cannot be used to compute QNM frequencies in 
such a case since there is an additional local minimum in the potential. Instead, we 
did an analytical study to show that an initial burst of radiation will settle into a 
constant value at late times. Such behavior have been observed for the I = 0 mode 
for the d Schwarzschild-de Sitter black holes. 

In addition to the WKB approach, we have used the unstable null geodesics ap¬ 
proach to compute u for massless scalar fields. First we observe that the radius of the 
unstable null geodesic Vnuii decreases with Q. The coordinate angular velocity, which 
is proportional to ur increases with Q. The Lyapunov exponent which is proportional 
to ojj does show similar behavior as was done using WKB approach. 

Finally, we have demonstrated that for the near-extremal black hole with Vf, ^ 
Tc the scalar field equation will have the Pdschl-Teller potential; one can get exact 
expressions for he QNM frequencies for this observation. Most of the observation we 
did with WKB approach, tally with the exact expressions obtained via this method. 

As future work, it may be interesting to compute QNM frequencies and study 
stability of the massive helds. We did not use WKB approach to compute QNM 
frequencies of this case since there may be extra local minima for the potential. One 
may use continued fraction method used by Konoplya and Zhidenko |4H] to compute 
massive QNM for the Schwarzschild black hole. 

Acknowledgments: The author wish to thank R. A. Konoplya for providing the 
Mathematica hie for the WKB approximation. 
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